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Abstract 

The bound-state method employed by Guo and Sheng (loc. cit.) is shown inadequate 
since only one of their solutions remains compatible, in the spin-symmetric low-mass 
regime, with the physical boundary conditions. We clarify the problem and construct 
a new, correct solution in the pseudospin-symmetric regime. 
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Recent letter PQ paid attention to the radial Dirac s— wave bound-state problem 
considered in the two different dynamical regimes. This problem has been reduced 
to the solution of the respective Schrodinger-like bound-state equations 

d 2 

- — F{r) + {M + E-C)[M-E + E{r)]F{r)=0, F(0) = F(oo) = 0, (1) 
(cf. eqs. (5) or (9) in [T]) and 

--^G(r) + (M-E + C)[M + E-A(r)]G(r) = 0, G(0) = G(oo) = 0, (2) 

(cf. eqs. (6) or (29) in [I]) where M denotes the mass, E is the bound-state energy, C 
is a parameter and the functions and A(r) represent certain phenomenological 
external potentials which have been chosen in the elementary Woods-Saxon form, 

E(r) = —, r , A(r) = r . (3) 

l+exp(^)' l+exp(^i) 

In the non-relativistic context, the solution of a very similar differential equation has 
been described in fair detail in ref. [2]. Unfortunately, the existence of certain specific 
physical features of the relativistic eqs. (JTJ and (J2J) forced one of us to imagine [3] 
that the constructions of ref. pQ should be re-examined. 

In a preparatory step the authors of ref. fTJ replaced the combinations of constants 
appearing in eqs. (JTJ), (J2J and © by the life-simplifying abbreviations 

= (M — E)(M + E — C)a 2 , fi 2 {ii) = (M + E){M - E + C)a 2 , 

v 2 {i) = (M + E - C) S a 2 , v 2 {ii) = — (M — E + C)A a 2 . 
This converts both eqs. and (j2J to the same eigenvalue problem, 
d 2 v 2 

- q2 ^2 y( r ) + ^ 2 y( r ) - x + exp [I. _ ^w Q i ^( r ) = °» y(o) = y(oo) = o. (4) 

The elementary change of variables 

r^x = 1/(1 + exp[(r — i2)/a]) 

reduces eq. to the Gauss's hypergeometric differential equation (compare eq. (64.3) 
in [2] with the same formulae preceding eqs. (9) and (29) in [I]). One arrives at the 
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exact wave functions sampled, say, by the unnumbered equation preceding eq. (64.10) 
in the non-relativistic construction Its relativistic analogue 

if>(x) = Nx^ [(1 - x) 5 /(//, 5) 2 F 1 (fi + 5, fi + 5 + 1, 1 + 25; 1 - x) 

+ (1 - x)- 5 f(fi, -5) 2 F 1 (fi - 5, fi - 5 + 1, 1 - 25; 1 - x)] (5) 

contains a certain ratio /(/i, 5) of products of pairs of T— functions given in full detail 
in eq. (13) of PQ. In this notation, the boundary condition in the origin becomes 
reduced to the single and exact transcendental equation 

* L ; =0- (6) 



1 + e~ R / a 

Usually [2], the purely numerical exact quantization condition (jUJ) is being replaced 
by its approximate asymptotic simplification 

V>(l-e- fl /°) • [l + 0(e- R/a )\ =0, i2/a»l. (7) 

The authors of ref. pQ did the same [cf. their four main quantization conditions 
(17), (22), (32) and (37)]. In this sense, the slightly misleading first line of their 
abstract (saying that their ". . . Dirac equation is solved exactly . . . " ) should rather 
read ". . . Dirac equation is solved approximatively . . . " . 

(i) In the first, spin-symmetric Dirac's dynamical regime the authors of ref. PQ 
had to separate the solution of the ordinary differential Dirac bound-state problem 
(JU) in two subcases. Incidentally, the first subcase [to be numbered or indexed as 
(i.a) in what follows] provides the only case where the result of ref. pQ is correct. 
We shall show that their construction fails to incorporate the boundary conditions 
properly in all the other cases. 

(i.a) Assuming that fi 2 ^ = ft 2 > and returning to eq. (9) of pQ in the first 
subcase one deduces that v%> = v 2 > and requires that the potential is sufficiently 
attractive, i.e., that v 2 —fi 2 = 7 2 > 0. In such a case, strictly speaking, the extraction 
of the exact spectrum from its implicit definition represents a straightforward, 
albeit purely numerical task. Unpleasant due to the well known slowness of the 
convergence of the infinite hypergeometric series 

_, , , . . ab a(a + 1)6(6 + 1) 2 
aFa a, b, c; z = 1 + — z + V - ' z 2 + ... 

c 2c(c + 1) 



in ip(x). Fortunately, in eq. (jHJ the value of the argument z remains exponentially 
small whenever R/a ^> 1. This means that in eq. (J7J) re- written in its fully explicit 
form 

~e~ R/a / (l + e- R/a ) 

= -/(//, -5) ■ aft (fi - 5, n - 6 + 1, 1 - 26; e- R ' a / (l + e"^)) 

the approximation 2 Fi (a, b, c; z) — > 1 will lead to very good precision. The authors of 
ref. |1| paralleled the non-relativistic exercise of ref. |2] and obtained the definitions 
(20) and (21) of the two components F n (r) and G n (r) of the related Dirac's bound- 
state wave functions. Their parallel implicit formula (19) for energies E n 

argT(2i7) — argT(/i + i7) — arctan(7//i) + jR/a = (n + l/2)7r (8) 

was simply based on the approximation (J7|). 

Let us note here that the next-order truncation 2^1 {a, b,c; z) — * 1 + abz/c clearly 
and manifestly confirms the expectations as it leads to an elementary explicit incor- 
poration of the next-order corrections via transition from eq. (jSJ) to another closed 
formula, 

7 R ( H 2 + l 2 \ ( R\ ( 1\ 

arg T(2i7) -arg Yiii+vy) -arctan — h7 — V \u + j+ —= exp = (n + - 7T. 

fi a \ 1 + 4y J V a J V 2/ 

(9) 

This improvement illustrates, persuasively enough, the exponential smallness of the 
higher-order corrections under the relativistic kinematics. 

(i.b) In an alternative range of parameters one admits that both the values of \i 
and v may be imaginary, i.e., /i 2 ^ = fi 2 < and vhs = v 2 < while [i 2 — v 2 = 5 2 > 
(cf. fP). It is fairly elementary to notice that at the very large r>lwe only have to 
deal with the significantly simplified constant-coefficient form of the original radial 
Dirac equation ((TJ), 

d 2 

- a 2 — F(r) + /4) F(r) = 0, r » 1. (10) 

This equation obviously possesses just the non- localized, asymptotically oscillatory 
solutions (remember that /1 is now purely imaginary). As a consequence, all the 
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"bound-state" solutions obtained in pQ remain non-localized and correspond in fact 
to the scattering states. In the other words, the implicit formula (24) of PQ for the 
energies E n is incorrect while the respective definitions (25) and (26) of the related 
Dirac's "bound-state" wave functions remain unnormalizable and, hence, irrelevant. 

An explanation of such a failure of the method returns us to the transition between 
formulae (11) and (12) in PQ where the latter formula has been erroneously declared 
to be the only solution compatible with the correct asymptotic bound-state boundary 
conditions. In fact, such an argument fails completely when \x becomes imaginary. 

(ii) Let us now move to the second, pseudospin-symmetric regime where the 
authors of ref. pQ claim the existence of the bound-state solutions of the Dirac bound- 
state problem (J2J) in the two separate subcases again [let us mark them as (ii.a) and 
(ii.b) in what follows]. We arrived at a different conclusion. 

(ii.a) In eq. (J2J) let us follow ref. pQ and assume, firstly, that /jfa* = fi 2 < while 
u fu) = ^ 2 > and v 2 — fi 2 = 7 2 > 0. This leads to the problem 



which parallels eq. IjlOj) . Hence, our conclusions remain the same. The implicit 
formula (34) for the bound-state energies E n in 1] is incorrect while the related 
wave functions (35) and (36) are not normalizable and represent merely a random 
sample of the scattering states. 

(ii.b) In the last subcase, equation (J2J) with the alternative choice of the real \i 
(with positive square = fi 2 > 0) is to be combined with the purely imaginary v 
such that z//U = v 2 < while /i 2 — v 2 = 5 2 > 0. Once we try to re- analyze such a 
model in present setting, we simply imagine that we deal with the Schrodinger-like 
equation 



where k 2 = M — E -\- C > and M + E > in the original notation. There are no 
bound states for the positive Ao > 0, of course. This observation is not incompatible 
with the observation made in ref. pQ where the wave functions (39) and (40) with 
the most elementary energy formula (38) represented just the single state lying at 
the very boundary \x = of the domain, anyhow. Nevertheless, a return to the latter 




G(r) + /4 t) G(r) = 



r > 1 



(11) 




(12) 
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/i = solution of eq. (J2J) "on the boundary" reveals that it does not need to exist at 
all. For example, we may verify that the insertion of the parameters \i = v = in 
the general solution (say, of the form similar to eq. (11) in [T]) leads to the exact 
wave function which is a constant. This constant must vanish due to the boundary 
conditions. 

In our final remark we have to emphasize that the situation in the pseudospin- 
symmetric regime is by far not so hopeless as it seems to be. As long as our potential 
functions are always asymptotically vanishing, our last equation (|12|) will possess (a 
finite number of) bound states whenever the effective potential well proves, paradox- 
ically, sufficiently strongly repulsive, — A(0) = +A < Aq Ulcal > < 0. 

The existence of the latter family of bound states was not considered in ref. 
at all. We repeat that in a search for a new relativistic Woods-Saxon model, one 
has to opt for a "paradoxical" choice of the repulsive barrier with A < (i.e., with 
A(r) > at all r). This choice, obviously, represents the only eligible pseudosp in- 
symmetric model where we would be allowed to construct bound states. The details 
of such a construction are entirely straightforward and may be left to the readers. 
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